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In this paper, the exact transient quantum transport of non-interacting nanostructures is investi¬ 
gated in the presence of initial system-lead correlations and initial lead-lead correlations for a device 
system coupled to general electronic leads. The exact master equation incorporating with initial 
correlations is derived through the extended quantum Langevin equation. The effects of the initial 
correlations are manifested through the time-dependent fluctuations contained explicitly in the exact 
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from the exact master equation. The resulting transient transport current can be expressed in terms 
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I. INTRODUCTION 

Quantum transport incorporating with initial corre¬ 
lations in nanostructures is a long-standing problem in 
mesoscopic physicsii In the past two decades, investiga¬ 
tions of quantum transport are mainly focused on steady- 
state phenomena^^— where initial correlations are not es¬ 
sential due to the memory loss effect. Recent experimen¬ 
tal developments allow one to measure transient quantum 
transport in different nano and quantum devices In 
the transient transport regime, initial correlations could 
induce different transport effects. In this paper, using the 
exact master equation approach j^^d^ we shall attempt to 
address the transient quantum transport incorporating 
with initial correlations. 

Conventional approaches for studying quantum trans¬ 
port include the scattering theoryi^iS the nonequilibrium 
Green function techniquejii^ and the master equation 
approachj ^h^^~ .l£ In the steady-state quantum transport 
regime, the famous Landauer-Biittiker formulcfc^ has been 
widely utilized to calculate successfully various trans¬ 
port properties in semiconductor nanostructures4i^ In 
the Landauer-Biittiker formula, the transport current is 
given in terms of a simple transmission coefficient ob¬ 
tained from the single-particle scattering matrix. How¬ 
ever, the scattering theory considers the reservoirs con¬ 
necting to the scattering region (the device system) to 
be always in equilibrium and electrons in the reservoir 
are always incoherent. Thus, the Landauer-Biittiker for¬ 
mula becomes invalid to transient quantum transport. 
The scattering theory method could be extended to deal 
with time-dependent transport phenomena, through the 
so-called the Floquet scattering theory;^ but it is only 
applicable to the case of the time-dependent quantum 
transport for systems driven by periodic time-dependent 
external helds. 

The nonequilibrium Green function technique based 


on Keldysh formalism^ has also been used extensively 
to investigate the steady-state quantum transport in 
mesoscopic systems, where the initial correlations were 
also ignoredfii^ In comparison with the scattering the¬ 
ory approach, the nonequilibrium Green function tech¬ 
nique provides a more microscopic picture to electron 
transport by formulating the transport current or more 
specifically the transmission coefficients in terms of the 
nonequilibrium electron Green functions of the device 
system in the nanostructures Wingreen et al. ex¬ 
tended Keldysh’s nonequilibrium Green function tech¬ 
nique to time-dependent quantum transport under time- 
dependent external bias and gate voltagesi^ But in the 
Keldysh formalism, nonequilibrium Green functions are 
dehned with the initial time to — oo, where the initial 
correlations are hardly taken into account, which makes 
the technique to be useful mostly in the nonequilibrium 
steady-state regime. 

The transient quantum transport was first proposed 
by Gini/3 under the so-called partition-free scheme. In 
this scheme, the whole system (the device system plus 
the leads together) is in thermal equilibrium up to time 
t = 0, and then one applies the external bias to let 
electrons flow. Thus, the device system and the leads 
are initially correlated. Stefanucci et. adapted 

nonequilibrium Green functions with the Kadanoff-Baym 
formalisnii to investigate the transient quantum trans¬ 
port with the partition-free scheme proposed by Cinii^ 
They obtained an analytic transient transport current 
in the wide-band limit. But in these works /^i^^ the 
transport solution is given with the nonequilibrium Green 
functions of the total system, rather than the Green func¬ 
tions for the device part of the nanostructureThus the 
advantage of simplifying the microscopic picture of quan¬ 
tum transport in terms of nonequilibrium Green func¬ 
tions of the device system may not be so obvious as in 
the Meir-Wingreen formulaJ^ 

On the other hand, the master equation approac h^h^^ 
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concerns the dynamic properties of the device system in 
terms of the time evolution of the reduced density ma¬ 
trix p{t) = TrE/3tot(t), where Trs is the trace over the 
environmental degrees of freedom. The dissipation and 
fluctuation dynamics of the device system induced by the 
reservoirs are fully manifested through the master equa¬ 
tion. Also, all the transport properties can be derived 
from the master equation. In principle, the master equa¬ 
tion for the quantum transport can be investigated in 
terms of the real-time diagrammatic expansion approach 
up to all the orders.— However, most of master equations 
are obtained by the perturbation theory up to the second 
order of the system-lead couplings, which is mainly appli¬ 
cable in the sequential tunneling regimej^ An interesting 
development of master equations in quantum transport 
systems is the hierarchical expansion of the equations of 
motion for the reduced density matrix,— which provided 
a systematical and also quite useful numerical calculation 
scheme to quantum transport. 

A few years ago, we derived the exact master equa¬ 
tion for non-interacting nano-deviceS )^^'^^ using the 
Feynman-Vernon influence functional approach^ in the 
coherent-state representation. The obtained exact mas¬ 
ter equation not only describes the quantum state dy¬ 
namics of the device system but also takes into account 
all the transient electronic transport properties. The 
transient transport current can be directly obtained from 
the exact master equation, which turns out to be ex¬ 
pressed by the nonequilibrium Green functions of the de¬ 
vice system. This unifies the master equation approach 
and the nonequilibrium Green function technique for 
quantum transport. However, the exact master equation 
given in Ref. (lail^ is derived in the partitioned scheme 
in which the system and the leads are initially uncorre¬ 
lated. The result transport current is consistent with the 
Meir-Wingreen formula.- This new theory has also been 
used to study quantum transport (including the transient 
transport) for various nanostructures recently.— >21 How¬ 
ever, realistically, it is possible and often unavoidable in 
experiments that the device system and the leads are 
initially correlated. Therefore, the transient transport 
theory based on master equation that takes the effect 
of initial correlations into account deserves of a further 
investigation. 

In this paper, we obtain the exact master equa¬ 
tion including the effect of initial correlations for non¬ 
interacting nanostructures through the extended quan¬ 
tum Langevin equation^ We find that the initial correla¬ 
tions only affect on the fluctuation dynamics of the device 
system, while the dissipation dynamics remains the same 
as in the case of initially uncorrelated states. The tran¬ 
sient transport current in the presence of initial system- 
lead and lead-lead correlations are also obtained directly 
from the exact master equation. This transient trans¬ 
port theory is suitable for arbitrary initial states between 
the system and its reservoirs, so that it naturally covers 
both the partitioned and partition-free schemes studied 
in previous works i^'^^’^^ Taking an experimentally realiz¬ 


able nano-fabrication system, a single-level quantum dot 
coupled to two one-dimensional tight-binding leads, as 
an specific example, we examine in detail the initial cor¬ 
relation effects in the transient transport current as well 
as in the electron occupation of the device system. 

In particular, we find that the electron occupation 
and transport current in the transient and also in the 
steady-state regimes can be different for partitioned and 
partition-free schemes when the system-lead couplings 
become strong, where localized bound states will ap¬ 
pear in the dot system. The quantum transport in the 
presence of localized bound states has also been studied 
in the literatures— In fact, Dhar and Sen consid¬ 
ered a wire connected to reservoirs that is modeled by 
a tight-binding noninteracting Hamiltonian in the parti¬ 
tioned scheme^i and they gave the steady-state solution 
of the density matrix and the current. Their results show 
that the memory effects induced by the localized bound 
states can be observed such that the density matrix of 
the system is initial-state dependent. Stefanucci used 
the Kadanoff-Baym formalism to formally study the lo¬ 
calized bound state effects in the quantum transport in 
the partition-free scheme— He found that the biased sys¬ 
tem with localized bound states does not evolve toward 
a stationary state. In this paper, we present the detailed 
time evolution of the density matrix and the transient 
transport current for both the partitioned and partition- 
free schemes. The initial-state dependance of the density 
matrix in the partitioned scheme is obvious in our formal¬ 
ism, as given in Eq. (HU. Meanwhile, the result of the 
system being unable to reach a stationary state found by 
Stefanucci^^ is also reproduced in our calculation. More¬ 
over, we find that this result is sensitive to the applied 
bias, as we will show in the paper. 

The rest of the paper is organized as following. In 
Sec. H, the transient quantum transport with initial cor¬ 
relations is formulated for non-interacting device system 
coupled to leads through the master equation approach. 
In Sec. HI, the transport dynamics with or without the 
initial correlations are explored in detail with a quan¬ 
tum dot device coupled to two leads modeled by one¬ 
dimensional tight-binding chains. The initial correlation 
effects in the electron occupation and in the transient 
transport current are explicitly calculated. Finally, a 
summary and discussion are given in Sec.IV. 


II. MASTER EQUATION WITH INITIAL 
CORRELATIONS 

To study transient quantum transport incorporating 
with initial correlations, we utilize the master equation 
approach. Consider a nanostructure consisting of a quan¬ 
tum device coupled with two leads (the source and drain), 
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which can be described by a Fano-Anderson Hamiltonian, 
H{t) =Hs{t) + HE{t) + HsE{t) 

~ 'y A 'y '^ £ak{t)c^^i.Cak 

ij ctk 

[^iakit)alcak + ^iak{^)^ak^'‘] ’ ( 1 ) 

iock 

where the electron-electron interactions are not consid¬ 
ered. Here, a\ (ai) and (cak) are creation (annihila¬ 
tion) operators of electrons in the device system and the 
lead a, respectively; (t) and eakit) are the correspond¬ 
ing energy levels, and Viak{t) is the tunneling amplitude 
between the orbital state i in the device system and the 
orbital state k in the lead a. These time-dependent pa¬ 
rameters in Eq. o can be controlled by external bias 
and gate voltages in experiments. 

Because the system and the leads are coupled only by 
the electron tunneling effect, and the leads are made by 
electrodes, the master equation describing the time evo¬ 
lution of the reduced density matrix p{t) = TrE[/5tot(t)] 
should have in general the following formi ^^d^i^^i^^ 

= - i[H'g{t), p{t)\ {7y (t)[2ojp(t)a| 

^3 

- aUjPi^) - P{t)4aj\ +%j{t)[a\p{t)aj 

- ajp{t)a\ + a\ajp{t) - p{t)aja3^ } 

= -i[Hs{t),p{t)\ [Ct{t) + C-{t)\p{t), 

cx. 

( 2 ) 


where the renormalized Hamiltonian Hg{t) = 
J2ij the coefficient eb(t) is the corre¬ 

sponding renormalized energy matrix of the device 
system, including the energy shift of each level and 
the lead-induced couplings between different levels. 
The time-dependent dissipation coefficients 'Jijit) and 
the fluctuation coefficients take into account all 

the back-actions between the device system and the 
reservoirs. The current superoperators of lead a, £+ (t) 
and C~ (t) give the transport current from lead a to the 
device system such that 


r /j-\ _ /dNa{t) ^ 

lUt)- e{ ) 

=eTr[£+(t)p(t)] =-eTr[£~(t)p(t)], (3) 


where Na{t) = J^k ^l,ki^)'^ak{t) is the particle number of 
the lead a. 


A. Without initial correlations 

When the whole system is initially partitioned, namely, 
Ptotito) = p{to) ® pE{to), in which the system can be 
in arbitrary initial state p(to) and the leads are initially 


at equilibrium pEi^o) = In this 

case, there is no system-lead and lead-lead correlations in 
the initial state. Then the time-dependent coefficients in 
Eq. ([5]) can be exactly derived using the Feynman-Vernon 
influence functional approach in fermion coherent state 

representation , with the results 

=^[u{t,to)u-^{t,to) -H.c.].. 

a 

= ^ XI( 4 b) 

Ct 

{t,to)v{t,t) -bH.c.]ij 

= - (4c) 


The current superoperators of lead a, and 

are explicitly given by 


^l{t)p{t) = - X {^«b(^)[alaiP(^) + alKOoi] 

■>-3 

+ Kaij{t)a\ajp{t) + (5a) 

^aii)p{t) =X + p{t)ajaW 

ij 

+ Kaij{t)ajp{t)al +II.C.} . (5b) 

In Eq. (jd]), u{t, to) and v{t, t) are related to the nonequi¬ 
librium Green functions of the device system in the 
nonequilibrium formalism.— These Green functions obey 
the following Kadanoff-Baym integro-differential equa¬ 
tions. 


dr 


t(r,to) + *eM(T,^o) + X / dT'ga{T,T')u{T',to) = Q, 

a -^*0 

(6a) 


d 

—v{T,t)+iev{T,t) + Y] / dT'ga{T,T')v{T',t) 

= X / dT'ga{T,T')u'<{t,T'), 

ry dto 


(6b) 


subject to the boundary conditions u(to,to) = 1 and 
v{to,t) = 0 with to < T < t. Here, the self-energy cor¬ 
relation functions from the lead to the device system, 
9 a{‘T,T') and ga{T,T'), are found to be 


9a,,(r,r') =X4"-'c(^)^/ofc(T (7a) 

k 

9aij{T,T') = X 4^*afc('r)E/„fc(r')/a(eafc)e"*47 

k 

(7b) 
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In Eq. (Hh), fa{eak) = + 1]“^ is the Fermi- 

Dirac distribution of lead a at initial time to- Solving the 
inhomogeneous equation Eq. (jBb) with the initial condi¬ 
tion v{tQ,t) = 0, we obtain 

/ dr' fdT"u{T,T')gaiT',T")u\t,T''). (8) 

In fact, Uij{t,to) = {{ai{t),aj{to)}) is the electron spec¬ 
tral Green function and v{T,t) is the electron correla¬ 
tion Green function of the device system in the non¬ 
equilibrium Green function technique.— The functions 
Ka{t) and Xa{t) in Eq. (HJ can be solved from Eq. ([6]), 

Kait)= / dTga(t,T)u{T,to)[u{t,tQ)]~^, (9a) 

Jto 

Xa{t) = / dT[ga{t,T)v{T,t) -ga(t,T)u\t,T)] 

Jto 

— Ka{t)v{t,t). (9b) 

The transient transport current flowing from lead a to 
the device system can be obtained directly from the cur¬ 
rent superoperators in the master equation and it 
is 

Ia{t) = - eTr[Aa(t) -|- Ka{t)(t) + H.C.] 

= -2eReTr/'dr[ga(t,T)p(^)(r,t) - gait,T)u\t,T)]. 
Jto 

( 10 ) 

In Eq. (flUl) , the single-particle correlation function of the 
device system pJ-\T,t) is given by 

= {a]{t)ai{T)) 

= [u{T,to)p^'^\U)u'{t,U) Jr v{T,t)\.., ( 11 ) 


our previously work^ we have used the extended quan¬ 
tum Langevin equation to correctly determine the time- 
dependent coefficients in the master equation in the ab¬ 
sence of initial correlations. Here, we shall use the same 
technique to determine the time-dependent coefficients 
in the master equation when the initial correlations are 
presented. The extended quantum Langevin equation 
for the device system operators can be derived from the 
Heisenberg equation of motion— 







-E 

Oij 


dTgaij{t,T)aj{T) 


*0 


i^H*afc(Ocafe(io)e (12) 


where the first term gives the evolution of the device 
system itself. The second term is the dissipation which 
arises from the coupling between the device system and 
the leads, and which is independent of initial state. As 
a result, the initial correlations cannot affect on dissipa¬ 
tion dynamics. The last term is the fluctuation induced 
by the reservoirs (the leads), which depends explicitly on 
the initial state of the whole system, so does initial corre¬ 
lations. The non-local time-correlation function g^ij (t, r) 
in Eq. (fT^ is given in Eq. ([7^), which characterizes the 
dissipation induced by the back-actions from the lead 
a to the device system. Since this quantum Langevin 
equation is derived exactly from Heisenberg equation of 
motion, it is valid for arbitrary initial state of the device 
system and the leads, including the case of the device 
system and the leads being initially correlated. 

Because of the linearity in the extended quantum 
Langevin equation (HU), its general solution can be ex¬ 
pressed as 


aiit) = ^ Uij{t,to)aj{to) + Fi{t), (13) 
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where p^^ (to) = (aj(to)az(io))) is the initial single¬ 
particle density matrix, and plj\T,t) = —iGfj(T,t) is 
indeed the lesser Green function in the nonequilibrium 
Green function technique, including a term induced by 
the initial occupation of the device system (i.e. the first 
term in Eq. CD)- Thus, the transient transport current 
obtained from the master equation has the exact same 
formula as that using the nonequilibrium Green function 
technique^ except for the above mentioned term origi¬ 
nated from the initial occupation of the device system 
which was not considered in Ref. 0. 

B. Including initial correlations 

When the device system and the leads are initially cor¬ 
related, i.e. ptotito) ^ p(to) (8> Psito), it would be chal¬ 
lenging to use the Feynman-Vernon influence functional 
approach to derive the master equation. Alternately, in 


where Uij (t, to) is the single-particle electron propagating 
function, and Ei (t) is an color-noise force induced by the 
electronic reservoirs (the leads). From Eq. (jT^ . one can 
find that Uij (t, to) and Fi(t) must obey the following 
equations, 

^Uij(t,to) +i[eU(t,to)]^^ 

+ E / dT[ga(t,T)U(T,to)].^ =0, (14a) 
a J*o 

d j-t 

—Ei (t) -f i ^ {t)Fj it)+Y^ / dTgo^ij (t, t)Fj (t) 

dt ^ Jto 

= -i ^ Viak (t)Cak {to)e~" -^‘0 ^ 

<yk 

(14b) 

subject to the initial conditions Uij (to, to) = 5ij and 
Fiito) = 0. It is easy to prove that the propagating 
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function Uij (t, to) is equal to the spectral Green function 
of Eq. (|nti): 

= ({a*(t),a](to)}) = (15) 

and the solution of Eq. (ITlb l is given by 

Fi{t)=-i^ fdT[uij{t,T)VjakiT)e~"-^*o 
akj 

(16) 

To determine the time-dependent coefficients in the 
master equation ©, we compute the equation of mo¬ 
tion of the single-particle density matrix of the device 
system, = (a](t)ai(t)) = Tr[a]aip(t)]. From the 

master equation (ED, the result is given by 

={p^^\t) [*e'(0 - lit)] 

- {[t£'it)+lit)]p^^Ht)},j+lijit)- (17) 

On the other hand, this equation can also be derived 
from the quantum Langevin equation. Explicitly, the 
single-particle correlation function of the device system 
calculated from the solution Eq. dH is given by 

p\fiT,t) = (a](t)a,(r)) 

= [uiT,to)p^'^\to)ut{t,to)]^.+vl{T,t). (18) 


Here, 


«ii(T 0 = t^ijiT,t)+[xiT,to)ut {t,to) + u{T,to)xHt,to)]^, 

(19)^ 


with 


pt 


XijitM) ={a]ito)F,{t)) = / dT[u{t,T)g^{T,to)]^J, 

a "tto 

(20a) 

t^^Ji^-,t) ={Fj{t)Fi{T)) 

a '^^0 to 

(20b) 


yOilJ 


(r,to) =iy^,ViakiT)e 


<^ck{ri)dTi, f 


a]ito)Cakito)), 


(20c) 


ti, (r, r')=Y.Il r... (r)e-‘« 


kk' 


X ^/c'fc' iT')e ^‘0 {cl,,,, (to)Cafe(to))- 

(20d) 


Physically, g^{T,to) characterizes the initial system-lead 
correlations, and g„(r, r') is a time-correlation function 
associated with the initial electron correlations in the 


leads. Compare Eq. dD with Eq. d , as one can see, 
the initial correlations modify the last term t;(r, t) in 
Eq. (ITTp by u'(r, t). Eq. (fT^ indeed gives the exact so¬ 
lution of the lesser Green function incorporating with 
initial correlations. 

From Eq. dUD, we can find 

^pij^(^) =[u{t,to)u-^{t,to)p^^'>{t)+R.c.],j 

- [u{t,to)u~^{t,to)v'{t,t) + H.c.]y 

+ j/^,it,t). (21) 

Now, by comparing Eq. (dD with Eq. EID , we can 
uniquely determine the time-dependent renormalized en¬ 
ergy e'ijit), the dissipation, and fluctuation coefficients 
7 y (t) and 7i_;(t) in the master equation incorporating 
with initial correlations. The results are given as follows. 


W = \ [^itM)u ^(t,to)-H.c.]^^. 

= Sij{t) — — ' [K,a{t) — kilit)]ij, (22a) 

OL 

lijit) = -^[uit,to)u-^it,to) +li.c.].^ 

= 2 ^ [^a(0 + k^ceit)]ij, (22b) 

O' 

Izjit) = - [u{t,tQ)u-^{t,to)v'{t,t)+R.c.]ij 

= + (22c) 

a 


From the above results, one can see that the renormalized 
energy and the dissipation coefficients are independent 
of the initial correlations and are identical to the results 
given in Eq. (jlDi) and (|41 d) for the decoupled initial state. 
The fluctuation coefficients also have the same form as 
that in the uncorrelated case, but the function v{t,to) 
which characterizes thermal fluctuation of the leads in 
the uncorrelated case is now modified by v'{t,to) which 
involves both the initial system-lead and the initial lead- 
lead correlations. In other words, initial correlations only 
contribute to the fluctuations dynamics of the device sys¬ 
tem. 

Correspondingly, Ka{t) in Eq. (1^ is the same as 
Eq. ([5^), but Xait) is modified by the initial correlations: 

Kit) = [ dT[gait,T)v'{T,t)-g',,{t,T)ut{t,T)] 

Jto 

- Kait)v'{t,t) + g^(t,tQ)u\t,to) (23) 
Thus, the transient transport current Ia{t) incorporating 
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with the initial correlations is given by 

lait) = - eTr[A^(t) + Ka(t)+ H.C.] 

= -2eReTr[/' dr[g„(i, (t, t) 

Jto 

- g'Jt,T)u'^{t,T)] + g^{t,toWit,to)]. (24) 

where g'^it, r) = J^a' Eq. ([Ml) contains all the 

information of transient transport physics incorporating 
with initial correlations, and is purely expressed in terms 
of nonequilibrium Green functions of the device system. 
It is easy to see that, if the whole system is in the parti¬ 
tioned scheme, then, g^{t, to) = 0 and g'^it, r) is reduced 
to ga{t,T) in Eq. dZb)- The single-particle correlation 
function Eq. (IT^ goes back to Eq. (HD, and Eq. (l24ll 
recovers the transient transport current Eq. HD for the 
partitioned scheme. 

In conclusion, the master equation ([2]) with the time- 
dependent coefficients of Eq. (1^ can be used to describe 
the non-Markovian dynamics of nano-device systems cou¬ 
pled to leads involved initial system-lead and lead-lead 
correlations. It should be pointed out that if the leads 
are made by superconductors where there may be initial 
pairing correlations. Then, the form of Eq. m may need 
to be modified further. We leave this problem for fur¬ 
ther investigation. Nevertheless, the master equation ([5]) 
is sufficient for the description of the transient quantum 
transport of nanostructures with initial correlations. The 
memory effects, including the initial-state dependence, 
are fully embedded into these time-dependent dissipation 
and fluctuation coefficients in our master equation. 

C. Relation to the nonequilibrium Green function 
of the total system 

In this subsection, as a consistent check, we prove 
that Eq. (IT^ can be reexpressed to the result obtained 


by Stefanucci and Almbladh in terms of the nonequilib¬ 
rium Green functions of the total system (i.e. Eq.(4) 

in Ref. HI). As we have shown, Uij(t,to) is the spec¬ 
tral Green function, Uij{t,to) = {{ai{t), Oj(to)}), the re¬ 
tarded and advanced Green functions of the device sys¬ 
tem can be expressed by Uij(t,to)'. 

Gfj{t,to) = -i9{t - to){{ai{t),a'l{tQ)}) = 

(25) 

Gfjito^t') =ie(t' -to){{ar{to),a]{t')}) =m^(f',to)- 

(26) 

Here, we have used the fact that the step function 9{t — 
to) — 9{t'—to) = I because t and t' are always larger than 
initial time to- The initial single-particle density matrix 
is associated with the lesser Green function of the device 
system: p[j\to) = -iG<{to,to) = (a](to)ai(io))- On 
the other hand, the nonequilibrium system-lead Green 
functions can be expressed by Uij(T,T'): 

G^ak(.t,to) = -i9{t - to){{a^{t),clf^{to)}) 

= -^[ dr'Mii(f,r')Via/c(r')e“*'^*o <^a.kiri)dTi 
I *^^0 

(27a) 

Gak,jito,t') = i9{t' -to){{cakito),a]it')}) 

=-^ r C. (r')<-(t, r'). 

l ''*0 

(27b) 

Thus, it is straightforward to find that 


=i[u{t,to)p^'^'’ {to)u'> {t' ,to)]ij + i[u{t,to)xHt' Go)]tj + i[x{t,to)u\t' ,to)]^J +iu,^{t,t') 
= ^ Gf,it, to)G<^, (to, to)Gf,^ [to, i') + E E to)G<^,ito, to)GXj ito,t') 


I'f 


i' ctk 


+E E Gi'^ak ^o)Gak,i' (^0 ) to)Gftj {to, ^0 + E E Gf'^ak ^o)Gak,a'k' (^ 0 , to)Ga'k',j (^ 0 , t') 


i' a.k ctk Oi'k' 

=[G^it, to)G< (to, to)G^ito,t%j = G<{t, t’), 


(28) 


r 


with 


G^,akito,to) = i{clkito)ai'ito)), 
Gak,i'ito,to) = i{al,ito)cckito)), 


and 

Gak,a'k'ito,to) = i{cl,k'ito)cak{to))- (30) 

This shows that the single-electron correlation function 
obtained by the extended quantum Langevin equation 
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gives the same results as the lesser Green function ob¬ 
tained by Stefanucci in terms of the nonequilibrium 
Green function of the total system.— 

III. ILLUSTRATION OF INITIAL 
CORRELATION EFFECTS 

In this section, we utilize the theory developed in the 
previous section to study the effect of initial correlations 
on transport dynamics. To be specific, we consider an ex¬ 
perimentally realizable nano-fabrication system, a single- 
level quantum dot coupled to the source and the drain 
which are modeled by two one-dimensional tight-binding 
leads. The Hamiltonian of the whole system is given by 

H{t) = £ca)a - ^(Aaio'l'cai -I- Xaicha) 

OL 

M 

+ ^ ^ ^ T ^a{t)\c]^n^an 

Q. n—1 

M-l 

~ ^ ^ ^ ' (Aac]j„CQ,„+i + (31) 

ol n—1 

where a (a^) is the annihilation (creation) operator of the 
single-level dot with the energy level £c, and Can {elm) 
the annihilation (creation) operator of lead a at site n. 
All the sites in lead a have an equal on-site energy ta- 
Ua {t) is the time-dependent bias voltage applied on lead 
a to shift the on-site energy. The second term in Eq. dSID 
describes the coupling between the quantum dot and the 
first site of the lead a with the coupling strength Aqi. 
The last term characterizes the electron tunneling be¬ 
tween two consecutive sites in the lead a with tunneling 
amplitude Aq,, and Af is the total number of sites. 

Using the Fourier transformation, the electron opera¬ 
tor at site n can be expressed in the k-space. 

Can = V^e““‘^^sin(nA;)cQ,fc, (32) 

k 

where fc = ^, / = 1 • • • A/", corresponding to the Bloch 
modes, and (j) = arg(AQ,). Then, the Hamiltonian (IHTl) 
becomes 

H{t) =£00^ a -I- £ak{t)c^akCak 

Oik 

+ (33) 

Oik 

where 

eak{^) — T Ua{i) j £ak — 2|Aq,| COS fc, (34a) 



The Hamiltonian (ESI) has the same structure as the 
Hamiltonian m in Sec. H. In the following, we consider 


two different initial states that are discussed the most 
in the literatures. One is partition-free scheme in which 
the whole system is in equilibrium before the external 
bias is switched on. The other one, which does not have 
any initial correlation, is the partitioned scheme in which 
the initial state of the dot system is uncorrelated to the 
leads before the tunneling couplings are turned on. The 
dot can be in any arbitrary initial state p(to) and the 
leads are initially at separated equilibrium state. Both 
of these two schemes can be realized through different ex¬ 
perimental setup. By comparing the transient transport 
dynamics for these two different initial schemes, one can 
see under what circumstances the initial correlations will 
affect on quantum transport in the transient regime as 
well as in the steady-state limit. 


A. Partition-free scheme 


In the partition-free scheme, the whole system is in 
equilibrium before the external bias voltage Ua{t) is 
turned on. The applied bias voltage is set to be uni¬ 
form on each lead such that Ua{t) = Ua&{t — to), so 
H{t < to) = H is time-independent. The density matrix 
of the whole system is given by 

Ptotito) = (35) 

where H and N are respectively the total Hamiltonian 
and the total particle number operator at initial time to- 
The whole system is initially at the temperature P = 
l/ksT with the chemical potential fj,. When t > to, we 
apply a uniform bias voltage on each lead. The whole sys¬ 
tem then suddenly change into a non-equilibrium state. 
When we take the site number Af to be infinity, the non¬ 
local time correlation function ga{t,T) can be expressed 
as 



Here the spectral density can be explicitly calculated 
from the tight-binding model with the result: 


ra(e) = 


?7a\/4|Aap - (e - Ca)^ if |e - Cal < 2|Ao 
0 otherwise. 


(37) 


and rja is the coupling ratio |Aai|/|Aa| of lead a. 

To calculate the initial system-lead and lead-lead cor¬ 
relations, {a'^{to)cak{to)) and {Ca,{to)cak{to)) , we diag¬ 
onalize the total Hamiltonian of Eq. (l33l) without the 
external bias through the following transformation 

«(io) = ^ VakG^{eak)bak, (38a) 

Oik 

Cak{to) =bak + '^ G^iCa'k') - 

e-a'k' - eak + 10 

Ol' k' 

(38b) 




with the retarded (advanced) self-energy, 


Thus, without applying external bias Eq. (1331) can be 
written as 


= X! ^akblk^ak- (39) 

Oik 

In Eq. (|38]) . 5 —>• 0+ and G^{eak) is the retarded Green 
function of the device system in the energy domain. The 
retarded (advanced) Green function is defined as 


S^’^(eafc) = 




Oi'k' 


^oik ^a.'k' ^ 


(41) 


G«’^(e„fc) = 


Ca/c Sc E^’^(eafc) 


(40) 


The initial system-lead and lead-lead correlations in the 
partition-free scheme are then given by 


{a\t^)cck{h)) =K:,G^(e„fe)/(e„fe) + ^^ 

Sa'k' - Sak + 10 

OL k 


I Vh'fe' \^G^{ea'k')G^ {Sa'k' )/ {Sa'k'), 


^a'k' (j'O^Oak (^o)) — bctk,a'k^ f (^afc) 




+ G {Sa'k')f{Sa'k') 


SoL'k' Sctk “t“ 


-|- G^ {Sak) f (Sak) 


Vc,,k'V* 


Oik 


^Oik ^ot'k' “t“ iS 


+ ^ \Va^ki\^G^{€aiki)G^{eaiki)f{Saikl) 




Oi'k' 




Oik 


aiki 


^ctifci ^oi'k' “1“ ^Oi'k' H“ 


(42a) 


(42b) 


r 


Using the initial correlations obtained above, we can cal- the non-local time system-lead and lead-lead correlation 
culate the initial density matrix in the dot, as well as functions. 


J 


p(i)(to) =(at(to)a(to)) = E |14fcpG«(e„fc)G^(e„fc)/(e„fc) = f 

gl{rM)=^ J |i[ReG^(e)]r„(e)f(e)e-H^+u<>)(--*»)+i J ^A(e)f(e) [?! ^ 

9lit,T)= J ^rc.ie)f{e)e-^^^+^Ait-r) 

+ X!/ —A{o)fie) [V J + J 

Ct' 

+ J 2 j ^[ReG^(e)]r„.(e)f(e)e‘^^+^“'^^"“*“^[iP J 


+ J2j ^[ReG^(e)]r„(e)f(e)e-'(^+U“)(‘-‘o)[p j ^^^^pi(d+u„,)(r-to 


~u)] 


“iE/ ^A(e)ra(e)r„,(e)/(e)e-*( 




(43a) 

(43b) 


(43c) 


where V represents the principal-value integral, and 
/(e) = -I- 1) is the Fermi distribution of the 

total system after diagonalized the total Hamiltonian. 
Thus, /(e) involves all the initial correlations in the 


partition-free scheme. 

To simplify the calculation, we made the band center 
of the source and drain match each other = sn = eo) 
and take Al = A/j = Aq to be the unit of energy. Also, it 
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should be pointed out that the diagonalization of Eq. (l38l) 
is performed in the absence of localized states of the sys¬ 
tem. The localized bound states can be added back via 
the spectral function A{e) = — 2ImG^(e)4 Without ap¬ 
plying the external bias, the spectral function ^(e) is 
given by 

A{e) =2n ^ ZjSie - ej)0(ry^ - ? 7 ±) 

2ImS--'-ne)e(4|Aop-(£-eo)^) 

where rj^ = Vl Vr- Eq. (ITil) . the first term char¬ 
acterizes the localized bonnd state with energy ej lying 
outside the energy band when the total coupling ratio 
> ^ 7 ^) where 77 ^ = 2 =F is the critical coupling ra¬ 
tio, and A = Ec — eo- As long as the energy bands of the 
two leads overlap, there are at most two localized bound 
states. The amplitude and the frequency of the localized 
bound state are given by 


partitioned schemes, we consider an initial empty dot in 
the partitioned scheme. Thus, the initial correlations for 
the partitioned scheme are 


{a'^ (to)a{to)) = 0 , 

(47a) 

{a\tQ)Cakito) = 0 , 

(47b) 

/ k' (Po')^ak(do')} — dak,a'k'fa{,^ak') ■ 

(47c) 


In this case, the non-local time system-lead correlation 
function ^^(t, to) = 0 , and the non-local time lead-lead 
correlation function is given by, 

9l{t.r) = + ( 43 ) 

which is indeed just the first term in Eq. (I43b b as an 
incoherent thermal effect from lead a. 


C. Transient transport dynamics 




1 ( 77 " - 2)v^ 4(7?^ - l)|AoP + 

2 (772-l)v^4(772-l)|AoP-f A2 


e± = eo 


{rf - 2)A rfsjAijf - l)|AoP -I- A^ 

2(772 - 1) 2(772 - 1) 

(45b) 


The effect of localized bound states are manifested in 
the second term of the system-lead and lead-lead corre¬ 
lation functions Eq. (l4^ i and (l4^i. respectively. Thus, 
the effect of initial correlations will be maintained in the 
steady-state limit. The second term in Eq. is the 
contribution from the regime of the continuous band, 
which causes electron dissipation in the dot system. The 
retarded self-energy in Eq. is 

- eo) - ^v'dlAoP - (e - eo)^]- (46) 

The above calculations specify all quantities related to 
initial correlations. 


B. Partitioned scheme 

For the partitioned scheme, the dot and the leads are 
initially uncorrelated, and the leads are initially at equi¬ 
librium state PE{to) = Te“AleAfter to 
one can turn on the tunneling couplings between the dot 
and the leads to let the system evolve. In comparison 
with the partition-free scheme, we also shift each energy 
level in lead a by Ua to preserve the charge neutrality, 
i.e. Cak Cak + Ua, and let Pl = Pr = /3.— Then the 
non-local time correlation function ga{t, t) is the same as 
Eq. (l36l) with the same spectral density of Eq. (l37l) . The 
spectral Green function u(r, tp) has the same solution 
as the one in the partition-free scheme. To demonstrate 
the initial-state differences between the partition-free and 


The purpose of this subsection is to inspect how the 
localized bound states and the applied bias affect the 
transient transport dynamics for the partition-free and 
partitioned schemes. The transient transport dynamics 
can be rationally described by the density matrix (i.e. 
electron occupation in the single-level quantum dot) and 
the transient transport current, 

=u(t,to)p(i)(to)M*(t,to) + w'(t,t), (49a) 

Ia{t) =-2eRe[ f dT[ga{t,T)p^^^ {T,t) 

Jto 

-gait,T)u*{t,T)]+g^{t,to)u*it,to)]. (49b) 

For the partition-free system, initially, the whole system 
is in equilibrium with temperature ksT, and chemical 
potential p. The source and the drain of the partitioned 
system are prepared at the same temperature {kRTR = 
ksTR = ksT). The systems are said to be unbiased 
under the condition eVso = Pl — Pr = Ul — Ur = 0. 

The time evolution of the absolute value of the spec¬ 
tral Green function |it(t)| = |u(t,to = 0)| which describes 
electron dissipation in the dot system is the same for 
both the partition-free and partitioned schemes in the 
unbiased and biased cases, as shown in Fig. [T] For the 
unbiased case, the spectral Green function purely decays 
to zero when there is no localized bound state, namely, 
the total coupling ratio rp < 2 —When one localized 

bound state appears, i.e. 2 —< 77 ^ < 2 -|--|^, |M(f)| os¬ 
cillates with time and then approaches to a non-zero con¬ 
stant value in the steady-state limit. This phenomenon 
becomes stronger when two localized bound states occur 
for 77 ^ > 2 -I- 1 ^. One can find that |M(t)| will oscillate 
in time forever when two localized bound states occur 
simultaneously. 

When the system is biased, \u{t) \ decays slower in com¬ 
parison with the unbiased system if the localized bound 
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state has not appeared. But when one localized bound 
state occurs, |M(t)| decays much faster and eventually ap¬ 
proaches to zero. In other words, the effect of the local¬ 
ized bound state can be manipulated by the bias voltage. 
When two localized bound states appear, the situation is 
dramatically changed, |u(t)| not only decays much faster 
in comparison with the unbiased case, but also the oscil¬ 
lation in the unbiased case is washed out in the steady- 
state limit, and |u(t)| approaches to a constant value. 
Numerically, we found that the amplitude of one of the 
localized bound states will be dramatically reduced when 
the bias voltage is applied. In other words, the applied 
bias can suppress the effect of one of the localized bound 
states. 

The above different dissipation dynamics in the bi¬ 
ased and unbiased cases will deeply change the transient 
and steady-state electron occupation in the dot and also 
the transport current. Fig. [5] shows the electron occu¬ 
pation in the dot and the transient transport current 
lL{t) = lR{i) in the unbiased case for the partitioned 
and partition-free schemes. The partition-free system is 
initially at equilibrium so that the dot contains a frac¬ 
tion of an electron, while the dot is initially empty in the 
partitioned scheme. One can see that the effect of the ini¬ 
tial correlations vanish in the long-time limit when there 
is no localized bound state. The steady-state electron 
occupation is the same for the partition-free and par¬ 
titioned schemes. However, when one localized bound 
state appears, the electron occupation in the dot oscil¬ 
lates in time for both the partitioned and partition-free 
schemes, and approaches to different steady-state val¬ 
ues. This indicates that the initial correlations can af¬ 
fect on the density matrix of the dot both in the tran¬ 
sient and steady-state regimes due to the existence of 
localized bound state. When two localized bound states 
occur, it will generate an new oscillation with the fre¬ 
quency being the energy difference of the two localized 
bound states for the electron occupation and this oscilla¬ 
tion will be maintained even in the steady state^^ where 
the initial-correlation dependence becomes more signifi¬ 
cant, as shown in Fig. He) for the electron occupation. 
The corresponding transient transport current for the 
partition-free and partitioned schemes approach to the 
same value in a every short time scale regardless whether 
the localized bound states exist or not. The steady-state 
current approaches to zero if there is no localized bound 
state. The current will also oscillate slightly when one lo¬ 
calized bound state occurs and approaches to zero in the 
steady state. When both localized bound states appear, 
the current will keep oscillations around zero forever 
However, the initial-correlations dependence in the trans¬ 
port current is not as significant as in the electron oc¬ 
cupation in both the transient and in the steady-state 
regimes, and even can be ignored in the steady-state 
limit, as shown in Fig. H 

The time evolution of the electron occupation in the 
dot and the transient transport current for both the par¬ 
titioned and the partition-free schemes for the biased case 


are shown in Fig. [3] Compare Fig. [2] with Fig.jSj one can 
find that the applied bias restrains most of the oscilla¬ 
tion behavior in the electron occupation as well as in 
the transport current, except for the very beginning of 
the transient regime. Also, regardless of the existence of 
localized bound states, the electron occupation and also 
the transport current all approach to a steady-state value 
other than zero due to the non-zero bias. In other words, 
the localized bound state has a less effect on the elec¬ 
tron occupation and the transport current when a bias 
is applied. This is because, as we have pointed out in 
the discussion of FigHc), the applied bias can reduce 
significantly the amplitude of one of the localized bound 
states, which suppress the oscillation of the transport 
electrons between two localized bound states. However, 
the remaining localized bound state will result in a small 
different steady-state values for partition-free and par¬ 
titioned schemes only in the electron occupation. The 
corresponding transient current flow through the left and 
right leads are quite different when a bias voltage is ap¬ 
plied. In particular, the transient transport current in 
the right lead is positive in the beginning for the par¬ 
titioned scheme because the dot is initially empty, and 
it approaches to a negative steady-state value in both 
schemes. But the steady-state current is almost inde¬ 
pendent of the initial correlations as shown in the inset 
graphs in Fig. H These result shows that the initial cor¬ 
relation effects are not so significant in the steady-state 
transport current in comparison with the electron occu¬ 
pation. 

In the end, as a self-consistent check, we consider the 
case Ul = Ur = 0. For the partition-free scheme, the 
system should stay in equilibrium. Our result is given in 
Fig.H which is expected. For the partitioned scheme, be¬ 
cause of the appearance of two localized bound states, the 
electron occupation in the dot will keep oscillation around 
a steady-state value, and the system would never reach 
to equilibrium with the leads. Also, the current con¬ 
tinuously oscillates with time around zero value. These 
results agree with the results obtained in Ref. and 
Ref. [ 2 ^, in which one consider the effect of localized 
bound state in quantum transport for partition-free and 
partitioned schemes, respectively. In fact, these results 
also agree with the fact that Anderson pointed out in 
Anderson localization^^ namely, the system cannot ap¬ 
proach to equilibrium when localized bound states occur. 


IV. SUMMARY AND DISCUSSION 

In summary, we investigate the transport dynam¬ 
ics of nanostructured devices in the presence of initial 
system-lead and lead-lead correlations. By deriving the 
exact master equation through the extended quantum 
Langvein equation, the effect of initial correlations is ex¬ 
plicitly built into the fluctuation coefficients in the master 
equation. The transient transport current incorporating 
with initial correlations is deduced from the master equa- 
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TlL=r|R=0.5|?.ol nL=ilR=l^ol TlL=nR=1-5|>-ol 



FIG. 1: The absolute value of time-dependent spectral Green function |u(t)| at different coupling ratios (a) i]l = i1r = 0.5|Ao|, 
(b) riL = rjn = 1.0|Ao|, (c) rjL = rjR = 1.5|Ao| with zero bias eVso = fJ-L — I^r = Ul — Ur = 0 (magenta solid line) and a finite 
bias eVsD = Hl — {J^r = Ul — Ur = 3|Ao| (navy blue dash line). The energy level of the quantum dot Ec = 3|Ao|, and the band 
center of the two leads eo = 2.5|Ao|. For the unbiased case, hl = fJ-R = 2.5|Ao|, and Ul = Ur = 1.5|Ao|. For the biased case, 
HL — 4|Ao|, fiR = |Ao|, Ul = 3|Ao|, and Ur — 0. We take to = 0 


r|L=riR=0.5|>,o| 'lL=nR=l^ol til=t1r= 1 .5|>,ol 



|>^o|t |A,o|t |A,o|t 

FIG. 2: The transient electron occupation of the dot and the transient transport current for the unbiased case at different 
coupling ratios (a) rjL = rjR — 0.5|Ao|, (b) tjl = rjR = |Ao|, (c) rjL = rjR = 1.5|Ao| for the partition-free (blue-dash line) 
and partitioned (red solid line) schemes. The energy level of the quantum dot Ec = 3|Ao|, the band center of the two leads 
£0 = 2.5|Ao|. For the partitioned scheme the leads are prepared at = Mh = 2.5|Ao|, and ksTL = ksTR = 3|Ao|. For 
the partition-free scheme, the system is initially at equilibrium with jj, = |Ao| and ksT = 3|Ao|. The applied bias voltage 
Ul = Ur = 1.5|Ao| after to = 0. 


tion. Our transient transport theory based on master 
equation approach derived in this paper is suitable for 
arbitrary initial system-lead correlated state. 

In application, we consider two initial states that 
are commonly discussed in the literatures, one is the 
partition-free scheme in which the initial correlations is 


presented explicitly. The other one is the partitioned 
scheme in which there is no initial correlations. The two 
schemes show the same steady-state behavior in the small 
system-lead coupling regime. When the coupling be¬ 
tween the dot and the leads gets strong, localized bound 
states will appear, and these two schemes give some dif- 
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r|L=r|R=0.5|?tol lL=nR=l>‘ol nL=^R=1-Sl^tol 



FIG. 3: The transient electron occupation of the dot and the transient transport current for the biased case at different coupling 
ratios (a) r/L = Vr = 0.5|Ao|, (b) rjL — rjR — [Aol, (c) rjL = tir = 1.5|Ao| for the partition-free (blue-dash line) and partitioned 
(red solid line) schemes. The energy level of the quantum dot = 3|Ao|, the band center of the two leads eo = 2.5|Ao|. For the 
partitioned scheme the leads are prepared at = 4|Ao|, ^ir = |Ao|, and ksTR = ksTR = 3|Ao|. For the partition-free scheme, 
the system is initially at equilibrium with /r = |Ao| and ksT = 3|Ao|. The applied bias voltage Ul = 3|Ao| and Ur = 0 after 
to = 0 . 



|A,o|t |>.o|t 

FIG. 4: The transient electron occupation n{t) and the transient transport current Ia{t) at the coupling ratio rjL = rjR — 1.5|Ao| 
with no energy shift {Ul = Ur = 0) for the partition-free (blue dash line) and partitioned system (red solid line). Other 
parameters are fJ, — fJ,L = jJ^R = \Xo\ and ksT = ksTR = ksTR = |Ao| 
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ferent results both in the transient regime and in the 
steady-state limit. It shows that the initial correlation 
effects can not be ignored if the localized bound states 
exist. Besides, when the localized bound states occur, the 
device system can not approach to equilibrium with the 
leads, and the initial correlation effects become more sig¬ 
nificant. These initial correlation effects accompanied by 
the localized bound state could be restrained by applying 
a finite bias voltage. Indeed, we find that a finite bias can 
suppress the oscillation behavior in the electron occupa¬ 
tion and in the transport current induced by the simulta¬ 
neous existence of two localized bound states. Compare 
with electron occupation in the dot, the steady-state cur¬ 
rent is insensitive to the initial states, so the results are 
not so distinguishable between the partitioned and the 
partition-free schemes. Nevertheless, the transport cur¬ 
rents are quite different in the transient regime in these 
two schemes when the bias is applied. Although we only 
consider a noninteracting nanostructure as a model in 


the detailed exploration of the initial correlation effects, 
it is not difficult to extend to include electron-electron 
interactions in our theory because the coefficients of the 
master equation are expressed in terms of the nonequilib¬ 
rium Green functions which can at least perturbatively 
or numerically handle the effect of electron-electron inter¬ 
actions. We leave this problem for further investigation. 
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